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Abstract 

The potential flow of two-dimensional ideal incompressible fluid with a free surface 
is studied. Using the theory of conformal mappings and Hamiltonian formalism 
allows us to derive exact equations of surface evolution. Simple form of the equations 
helped to discover new integrals of motion. These integrals are connected with the 
analytical properties of conformal mapping and complex velocity. Simple form of 
the equations also makes the numerical simulations of the free surface evolution 
very straightforward. 

In the limit of almost flat surface the equations can be reduced to the Hopf 
equation. 
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1 Basic equations 



We study the potential flow of two-dimensional ideal incompressible fluid. The 
fluid occupies a half-infinite domain 

— oo < y < r](x, £), — oo < x < oo. (1.1) 
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The flow is potential, so that: 

u = V$, ^\ y=v{ ^ t) = ip(x,t). (1.2) 

Boundary conditions on the surface are standard: 



aT + 2 |v$l + ^ = P ' 

^ + ^ = $ y (1.3) 
at y — r)(x, t) and 

— = 0,7/ -»■ -oo, 
ay 

5$ 

- = 0,1400. (1.4) 
ax 

Here g is the gravity acceleration and P - constant pressure at the surface 
(Let P = 0). It is known ( jZakharov (1968)] ), that the shape of surface rj(x, t) 
and the potential on the surface ij){x,t) form pair of canonically conjugated 
variables obeying the Hamiltonian equations: 

dt~ Sip 1 dt ~ Srj' { ' ' 

Here H is Hamiltonian function (total energy of the fluid): 

U = T = -\ dx \V$\ 2 dy 

2* J —oo J — oo 

kinetic energy of the fluid. Along with the energy there are three more integrals 
of motion, amount of fluid: 

d f°° 

— J r](x,t)dx = 0, (1.6) 
and vertical and horizontal momenta: 



Q [oo rrj{x,t) 

dx I 

J — oo 

dx <j> x dy = 0. (1.7) 



/oo pri(x,t) 
dx / <f> y dy = 
-OO J — CO 

Q poo rrj(x,t) 



at 



oo J —oo 



The equations (11. 3p are functionally nonlinear and can be hardly studied. The 
most important known solutions were derived by Dirichlet in 1860. In these 
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solutions a shape of the surface is a quadric (ellipse, hyperbola and parabola). 
Dirichlet solutions are described in details in jLonguet-Higgins" (1976)] . 

Equations ( II. 5p minimize the action 

S = I Ldt (U 



L= ipr] t dx-H. (1.9) 

J — oo 

Starting from this point let us forget for a while about hydrodynamics, and 
consider more general case. Namely, let's think of % as some arbitrary func- 
tional of ip and 77. Then, let's introduce instead of rj(x, t) new complex variable 
z(u, t) in the following way: z(w, t) is complex function of complex variable 
w = u + iv. It is analytic in the lower half-plane of complex variable w. Its 
real and imaginary parts, given on the real axes (z(u, t) = x(u,t) + iy(u,t)) 
parametrically define y(x,t) = r](x,t). 

Another words, z(w,t) is the conformal mapping of the domain, bounded by 
the curve r](x, t) to the lower half-plane of w, 

w = u + iv, —00 < u < 00, —00 < v < 0. 



In the new variables Lagrangian takes the form 



/oo roo A 

ip[y t x u - x t y u ]du - f(y- H(x - u))du - H. 
-00 J — 00 

Here H is the Hilbert transformation, that provide the relation between real 
and imaginary parts of analytic function: 

y = H(x-u) = -P.V. r X ^~ U ' du'. (1.10) 



7T 

Motion equations should be found from the condition 5S = 0. When per- 
forming variation one should take into account that x and y are connected by 
relation ( ll.iop . and Hamiltonian function T-i can be considered as a functional, 
depending only on ip and y. Hence, one can put #f = 0. 

The condition 5S = leads to the following "implicit" equations of motion 

y t x u - x t y u = — (1.11) 
dip 

%p t x u - x t ip u + f = — r~ ( 1 - 12 ) 

by 
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VuA-yt^u-Hf = (1.13) 

From fll.l3f> one can find 

f = H{y t ij Ju -yM- (1-14) 
While ip satisfies the equation 



bl-i 

ip t x u - x t if) u - H{y u ip t - yM = — (1.15) 

by 

Equations ( II. lip and (I1.15P are resolved with respect to variational derivatives 
|^ and They can be written in the explicit Hamiltonian form: 



r , . ; sn 

ll n y + ^12^ = —, 
by 

-n£y=^. (1.16) 
Here the operators comprising the symplectic form f2y are: 



n u = -(^ + (^)), 

ft 12 = -a/ + iiy, fi£ = x' + y'H. 

As far as the equations f II . 1 6 j) are obtained directly from the variational prin- 
ciple, the symplectic form Qij is close and nondegenerated. It means, in par- 
ticulary that equations ( 11.1611 have no Casimirs - the constant of motion which 
do not depend on a choice of the Hamiltonian function %. 

This is remarkable that equations II. lip and fll.1511 in a general case can for 
an arbitrary % be resolved explicitly with respect to time-derivatives. To do 
this we introduce projective oprerators 

P ± = -(1 ± iH), P+P- = 0, P ±2 = P ± 

. Any complex- valued function (f)(u), — oo < u < oo can be presented as follows 

(j){u) = P + (f)(u) + P~(f)(u) 

. Here P ± <p(u) are analytic functions in the upper and lower half-planes. 
Mention first that equation (jl.lip is equivalent to the equation 

™ (1.17) 
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We denote the Jacobian J of the conformal mapping 

J=\zf. 

Dividing (I1.17P by J and applying projective operator 

P- = \(l + iH) 
one can get from (11.171) the following equation: 

z — illz', 

^2P-(-|) (1.18) 

For y and x one can get the equations: 

y=[x — y H) , 

y V y J Sip 

x = -{y' + x'H)- J °-^. (1.19) 

Excluding t/j t from (11.121) and (11.131) and using equation (II. lip one can obtain 
the folowing relation: 



y'f + x'Hf = -y'^ - ^ = -(Fz> - Fz>). (1.20) 

Here F = f + iHf, function being analytic in the lower half-plane. By dividing 
on J and applying the projective operator one can find: 



F -~7 P ^ + (L21) 

Note that expression (jl.2ip does not include time derivatives. Now one can 
express ip t from equations (I1.12p and (I1.13I) . using equations (II. lip and (ll.2ip . 
After simple calculations we end up with the following equation: 



*=-«'«4 + 7>^>» + ^ < L22 > 

Equations (II. 19ft and (11. 22ft can be written in the "implectic" form 
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dtp oy 

y=R^ (1.23) 



Here 



f?ll 


=-( 


R22 


= 0, 




1 


R\2 


= 7' 


R21 


= (x' 



A 1 

1 + T 



v'H)j,R 2l = fl+. (1.24) 

These formulae determine the Poisson structure on the functionals defined on 
the real functions y and ■0. Let a and /3 is a pair such functionals. Obviously 



{ 1 '} I °° { ^ a R ^ ^ a R ^ + ^ a R ^ }dx 

-00 Sip(x) 5ijj(x) 5ip(x) Sy(x) Sy(x) 5ij){x) 



2 Basic equations in the complex form 



It is convenient to accomplish equation (11.181) imposed to the complex poten- 
tial 



® = tp + iHip = 2P-il;, (2.25) 
$ = <0-iH^ = 2P + ip. 

Note that 
Hence 



[/ = 4P-{i(p-^ + P+^)}. (2.28) 
Now equation (11.181) can be presented in the complex form. 
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Now one can apply operator 2P to the equation ( ll.22p and express |^ and 
^ by the use of (I2.26j) and (12 . 27j) . We get closed equation for $. This equation 
can be transformed to the following simple form ( see Appendix [B] ) 

$ = iU& -B-V. (2.29) 

Here 

B = -4iP-{j(P-(®'^) ~ P + (& ji))} (2-30) 

and 

V = -Aip-{hp-z> 5 -lL - P+z'^)}. (2.31) 
J oz oz 

Equations f 1 1 . 1 8 j) and (I2.29|) compose close system of Hamiltonian equations 
written in the complex form. Poisson bracket in terms of (z, $) is discussed in 
Appendix IA1 



There is another form of complex equations. Following to the article Dyachenko (200f)] 
we introduce new variables: 

z 

V = i— = iR&. 

oz 

In terms of R and V equations (II. 18ft and (I2.29f) take a form: 

^ = i (UR>-RU>), 

^ = i(UV-R{B + V)'). 



It is important to stress that equations (11.181) . (12.291) and (12 .32|) are written 
for the functions which are analytical in the lower half-plane (Imu> < 0). 

We should stress once more that equations (I2.32j) are just another form of 
the general Hamiltonian equations (11.51) . Meanwhile, this particular form of 
the Hamiltonian equations with the very noncanonical Poisson brackets is 
in our opinion the most convenient for both analytic and numeric study of 
hydrodynamics with free surface. Finally we present implicit equations (II. lip . 
(I1.15P in the complex form. To do this we remember that 



x ■ 



2^ + ^)> 



V = — \ z — z), 
y 2V ' 



(2.32) 
(2.33) 
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One can see that equations (11.111) . fll . 15j) can be presented as follow 



P- 



A ,6H 

z t z u - z t z u + Ai—=- 
5$ 



§ t z u - Zt<$> u + P <& t z u - z t ^u + 4z 



0. 

1¥ 



o. 



(2.34) 
(2.35) 



3 Hydrodynamics with a free surface 



All the machinery described above was designed to study different versions of 
potential flows of ideal fluid with a free surface. We start syatematic study of 
these models from the simplest one - "free" incompressible fluid moving only 
under influence of force of inertia. In this case the Hamiltonian is just kinetic 
energy of the fluid 



1 r 00 j roo _ _ 

H = T=-- if>Hif/du = - + 

2, J —CO o J — CO 



(3.36) 



Now 



5U 



-Hip, (3.37) 



^ 1 ^ i 
5$ 4 5$ 4 

Let us consider first the implicit equation. Equations fll. lip and (I1.15P take 
the following form 

y t x u - x t y u = -Hip u , (3.39) 

ifi t x v - x t tp u - H(y u ip t - y t ip u ) = 0. (3.40) 
Equation (I3.39|) can be rewritten as follow 

z t z u - z t z u = -2iHip u . (3.41) 

This is the step to complex form of implicit equations. Full complex implicit 
equations are: 

z t z u - z t z u = $ u - $ u , (3.42) 
or after applying projector operator P~ 

P- {z t z u - z t z u ) = -$ u . (3.43) 

Second equation is 

<f> t z u - z t <$> u + P- (® t z u - zt$ u ) = 0. (3.44) 
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The "transport velocity" U is defined only by kinetic energy For all hydrody- 
namic model U and B are universal 



1 - \ ~ / $ — $ 



u =- 2P [j H ^)= iP {-w^w)> (3 - 45) 

B = p-( l pQi. (3.46) 



Thus "scalar" explicit equations read: 



y——{x' — y'H) ^jHipu, 



ip = i> u ^p + jH{ip u H^ u ). (3.47) 
The "complex" explicit equations has the standart form 



z — iUz', 

Q = iU&-B. (3.48) 

U and B are given by equations ( I3.45p . Finally the Dyachenko equations for 
this simplest case read 



R t = i(UR' - RU'), 

V t = i(UV' — RB'). (3.49) 
In R and V variables 



U = P-(RV + RV), (3.50) 



B = P-{VV). (3.51) 

One more usefull form of "free" hydrodynamic equations stems from (jl.2ip . 
Apparently 



F = -Lpr-W *) = -P-^'H^') = (3.52) 

z dtp z 

igP-d^) = £. (3.53) 
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From (11.12p and f l 1 . 1 3 j) we obtain important equation 



$2 

^z u - z^ u + -\ = 0. (3.54) 



This is the Bernoulli equation written in conformal variables. 



4 Dirichlet solutions 



"Free" hydrodynamic equations admit a remarkable class of solutions, found 
first by L. Dirichlet in 19 century. In These solutions potential is a quadratic 
function of coordinates, and the free surface is a curve of second order. It is 
interesing to find these solutions in conformal variables. Let us assume that 



$ = a ( t ) + M^ 



(4.55) 



a(t) and b(t) are real. Then 



$' = bzz', 



V>= ^($ + $) = a+ b -{z 2 + z 2 ) 



Equation (I3.42p reads now 

zz! — zz' = b(zz' — zz'). 
Equation (I3.54p takes the form 
b 



a + -(z 2 + z 2 
4 V 



z' + ^-(zz' -^z') + h -z 2 z' = 0. 



(4.56) 



(4.57) 



Using equation ( I4.56P we realize that the terms proportional to z' are cancelled. 
Then we cancel z and end up with the equation for the quadric 



b b 2 
a + -(z 2 + z 2 ) + —\z\ 2 
4 V ' 2 1 1 



0. 



(4.58) 



This equation describe the shape of the surfce. To find conditions for d, b we 
should solve the equation (14.561) . 

Suppose that b — b 2 > 0. Now quadric ( I4.58P is hyperbola 



^!-=i + — 

A 2 B 2 ' 

A 2 = J^, 
b-b 2 



B 2 



2 a 



b + b 2 



(4.59) 
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Equation (I4.56P can be rewritten as follow 

yx — xy = —b(xy' + x'y). 
By plugging ( 14. 59ft into (I4.60[) one get 



(4.60) 



A = -bA, B = bB, 

AB = X 2 = const, B 2 = — (4.61) 

A 2 

Then 

A 2 (b-b 2 ) = B 2 (b + b 2 ) 

or 

Kl-^) = ^ 2 (1 + ^) (4-62) 

From (I4.62p one can see that the case A 2 = B 2 = X 2 is marginal. If A 2 < A 2 , 
b < 0. in the opposite case A 2 > A 2 , b > 0. By assuming that 

■db A1 db 

b = ^TT = -^77, 
dA dA 



one get 

db 1 A A + X 4 



dA AA A - X A ' 
Thus the equation can be easily integrated 

FA _A 

|^4_ A 4|l/2 A ' 



(4.63) 



(4.64) 



F > - is constant of integration. Equation (14.641 can be integrated elliptical 
functions. In two limiting cases A « X and A > X it essentially simplifies. 

Let A\ t=0 = A . 
If A « A 



A = -FA 2 , A = —, q=—. (4.65) 

t + q A 

If A » A 



A — F, A = Ft + A . (4.66) 

Hence the "flat" hyperbola A < X becomes more flat, while the "sharp" 
hyperbola A > X becomes more sharp as t — > oo. 
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Let us introduce new variables c and such that 

A = c cos 0, B = c sin 9 



2 = _!£(i^ + I e -*>) ; 

2 v # ; ' 



(4.67) 



R = e farcsinh( W ) = ( w + v^^T) f 



(4.68) 



Then 



^ 29 



29 



R = — arcsinh{w)R = — (u> + Vw 2 + l)R, 

71 TV 



(4.69) 



R' — , =R, 

vr Vw 2 + 1 



(4.70) 



ic 1 , „ .a 1 
(Re ld + —e 



(4.71) 



[c + ic9)Re ie + (c-ic9)^e 

R 



— ■ — arcsinh(w)(Re ie + — (< 

TV R 



Now remember that 



AB = c 2 sin 9 cos = -c 2 sin 20. 

2 



(4.73) 



By differentiating (I4.73P one get 



c cos 20 • 

c= —9. 

sin 20 

Finally we got 



(4.74) 



7>0 112 I 

* = L___(ife* + e" ie ) - -arcszn/i( W )(i?e ie + -e" ie ) 
2 sm 29 R it R 



(4.75) 
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or 



sin 29 
As we know 



dotO _ it9 — - . 

z + — v!fl + larcsinn[w)z . 



(4.76) 



-c cos 9 + c9 sin (9 
ccos 



C0 



sin26> 



(4.77) 



Tj get the last result we need equation (14.631) . Now by plugging (14.721) , (14.761) , 
( I4.77P into (14.561) we see that this equation is automatically satisfied. Now we 
can find the transport velocity U. According to ( 13.481 ) 



U = - % ±. 
z' 



(4.78) 



By the use of (14711 and (14761) we get 

m9 



U 
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Vu> 2 + 1 



1 Re ie + ^e~ w , . ^ — - 
f— g + log (tu + Vu> 2 + 1) 



2sin2#f?e ie - \e 



(4.79) 



Note that 2 has only one branch point at w — i. It is analytic everywhere 
exept the cut along the imaginary axis (i, zoo). The same statement is correct 
for all other analytic in the lower half-plane functions 



R = -,V = bz.B = P-(VV). 



Apparently 



Note that 



R= ~9^ fTl Re*-^- 



1 



R\ w =i — 0, U\ l 
To find B we again can use equation ( I3.48p . 



(4.80) 
(4.81) 



5 This is the end... 



To apply the developed general theory to concrete physical system one should 
specify the Hamiltonian function H. We will study the following class of Hamil- 
tonian functions: 



T~L — Hi + H2 + ^3 + Hi- 



(5.82) 
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Here 



ory fOO fy POD 

-f {z 1 + z' -2){z-z)du = ~ / y{x'-l)du, (5.83) 

o J —CO 2i J — CO 



U 3 = -— r {z - zf{z' + = 9 - y 2 x'du, (5.84) 



16J-oo v 7 v 7 2 7- 



OC 



h i = a /°° I V^¥ - ^ * * j du = a J°° ( vV 2 + x ' 2 - x ') du - ( 5 - 85 ) 



Hamiltonian "Hi is the kinetic energy of the ideal incompressible fluid. It de- 
pends on $ only, not on z. The others, H2, H3 and I-L4, are different forms 
of potential energy depending on z only. Hamiltonians "H 3 and "H 4 are the 
components of the potential energy in the presence of gravity acceleration g 
and surface tension o. Hamiltonian "H 2 has more sophisticated physical in- 
terpretation (see Zubarev (2000)| ). For a < it is a potential energy of the 



dielectric fluid with ideally conducting free surface, placed into the electric 
field. For a > this is the potential energy of the same fluid containing inside 
the magnetic field. More detailed information is in the Apendix B. 

For the Hamiltonian function (I5.82p 



1 $' 1 $' 

U = zP"{- 7 - --} = P-{RV + RV), (5.86) 



£ = = p-(yV). (5.87) 

In the absence of a potential energy V = (see (12.3ip ). In the presence of a 
potential energy 



V = V 2 + V 3 + V 4 . (5.88) 

Here 

P 2 = aP-(RR-l), (5.89) 

Va = -ig(z-w), (5.90) 
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Vi = 2aP-(Q'Q - Q'Q), 



q = Vr. 



(5.91) 



In the absence of the surface tension motion equations (I2.32p take the following 
remarkably simple form: 



R t = i(UB! - RU'), 

V t = i{UV' - RP~(VV + aRR)') + g(R - 1). (5.92) 

Equations (I5.92p with (I5.86P are cubic with respect to the unknown functions 
R and V. It makes them very suitable for numerical simulation. 

Let us put now a = 0, and include into consideration the surface tension. The 
first equation of (15.921) is not changed, while the second takes the form: 

V t = i(UV' - RP-(VV)' + g(R - 1)) - 2aRp-(Q'Q - Q'Q)', (5.93) 
For function Q equations (I5.92p became quartic nonlinear: 

Qt = i{UQ - hj'Q), 

V t = i(UV - Q 2 P-(VV)') + g(Q 2 - 1) - 2aQ 2 p-(Q'Q - Q'Q)', (5.94) 
and U is equal to: 

U = P~{VQ 2 + VQ 2 ). 
Equations f 1 5 . 9 2 j) keep the same form for other boundary conditions. Just for 
different boundary conditions we need to replace Hilbert transformation by 
some other operator. For instance for periodic boundary conditions one has to 
use Schwartz transformation instead of Hilbert. Another important case, fluid 
of finite depth is discussed in |Dyachenko, Kuznetsov Zakharov| . 

Note that equations include now only derivatives of conformal mapping and 
complex velocity potential. As regards integrals of motion they acquire more 
complicated form. But if one restores complex velocity potential 

f V 

$ = —it —dw 
J R 

then the kinetic energy is equal to 

/oo 
Re($)Im(&)du 
-oo 

and momenta (given by (11. 7p ) are now equal to: 

Re($)Re( — )du, P x = / Re($)Im(—)du. 

-OO R J— oo R 



15 



6 New integrals of motion 



Equations ( I5.92p have some remarkable properties. Suppose that z(w) and 
have at w = A a logarithmic branch-point (ImA > 0): 



z = a log (w — A) + w, 

$ = Mog (w - A). (6.95) 
At this point 



2 = r + 1, 

w — A 

$' = — ^— (6.96) 
w — A 

and 

$' 6 1 
1/ = z — ~ i-. i? ~ -(w — A). 

z' a c 

One can see that both R and V are regular in the neighborhood of point A. 
the same statement is correct for B, P 2 and TV Note that it is not correct for 
V\. Indeed, near w = A 



(T 



— A 



Regularity of U and i? results in the following important consequence. Suppose 
that functions R and V are regular not only in the lower half-plane, but also 
in some domain above the real axis. And in this domain R has simple zeroes 
at the points 

Ai, A 2 , . . . A n . 

In the vicinity of each zero 



R~a n (w - A n ), 
and V ~V n . 



If R and V satisfy the equations (15.921) . then 



da n dV n 

~& = °' ~dT = g (6 - 97) 

In the other words, coefficients a n , which are nothing but logarithmic residues 
of z, are constant of motion of the equations (I5.92p . In the absence of gravity 
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logarithmic residues of $ are constant of motion either. In presence of gravity 
they are linear functions of time 



V n = V n0 + gt 

and 

b n = -ia n V n = -ia n (V n0 + gt). 

It is important to find the Poisson's bracket between new integrals of motion. 
As far as a n are completely defined by z, so that 

SOn _ SOn _ n 

one can conclude: 

{a n ,a m }. (6.98) 
Poisson's brackets 

{a n ,b m }, and {b n ,b m } 

are still unknown. 

So far we have discussed only simple zeroes of R. In fact, one can release this 
condition and consider zeroes of higher order. One can assume that in the 
vicinity of A 

U ~ m + (w — X)ui, 

R~(w-\) n R. (6.99) 



Plugging (I6.99P into the (15 .92 j) one can find 



R = i(UR! -RU' + nmR), 

\ = -u . (6.100) 

This assumption (16.991) is structurally stable. One can assume again that V, 
U, B and V are regular at w — A. For 

Vn = V\ w=Xn 

one can obtain the following formula: 



f = 9 . (6.101) 
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In this case both z' and $' have a pole of n-th order at the point w = A: 



, a- n a_i 

2 = 7 rT~ + ^ 7 77 + • • • ) 

(w - A) n (w — A) 

$ , = 7 -^ + --- + 7 -^T + ---- (6-102) 
(u> — A) n (u> — A) 

Plugging (I6.102p into the equations (I3.48p . one can obtain: 



d d 

— a_i = 0, — 6_i = a_i5i (6.103) 

at at 

Moreover, (16. 10 If) reads: 

d 0_n 



(it a_ 



^. (6.104) 



it 



Formulae (16. 1Q3[) can be written as follow: 



4- i z'dw = 0, 4- f ®' dw = gf z ' dw - (6.105) 

at Jr dt Jr Jr 

In ( I6.105P one integrates along small contour around A. 

Formulae (I6.10ip and (16. 105[) do not include the parameter a. One can think 
that they could be extended to a more broad class of the Hamiltonian func- 
tions. This is actually not true. Only a very special Hamilton functions keep 
U, B and V analytic at the points where z' and $' have poles. In a general case 
singularities in z' and $' generate singularities in U, B and V, and formulae 
flpUOTD . flBTTUBl . flBTTUD and flBTTUoD are violated. 

This phenomenon can be traced for the case when we include into considera- 
tion surface tension. If R has a simple zero 



R~(w — A), (6.106) 
the expression for surface tension term in (15.941) is 



-2aRP~(Q'Q-Q'Q)' ~ r (6.107) 

(w — A) 2 

It means that analyticity of V is immediately violated. The situation can be 
fixed if Q = y/R is analytic functions, having zero of any arbitrary integer 
order, one can see that the assumption 



18 



Q = (w - \) n Q 
is compatible with equations (I5.94p . 
Expression (I6.108[) for Q implies that 



(6.108) 



R = (w- \) 2n R. (6.109) 
Thus R has zero of an even order. 

In conclusion of this chapter we should stress that even in the absence of 
gravity integrals (16.1011) and (I6.1Q5[) cannot be interpreted as Casimirs for 
some degenerated Poisson's bracket. Just due to the fact of existence of the 
variational principle and the symplectic structure, the Poisson's bracket in 
this case is not degenerated, and no Casimirs do exist. 



7 Connection to the LGE equation 



In the articles Dyachenko fc Zakharov (1996)| , jZakharov fc Dyachenko (1996)] 



we have found that in certain situations the dynamics of a free-surface fluid 
can be described by the Laplace growth Equation (LGE), which is known 
in hydrodynamics since 1945 ( |Polubarinova-Kochina (1945)] , |Galin (1945)1 ). 
We show now that the LGE appears in a natural way in the framework of our 
formalism. 

Let us suppose that the fluid flow has two scales - global large scale with the 
complex potential \l/o; and small scale with the potential 5i/j. Another words: 

y = ^ + 6tp. (7.110) 

Let us denote 

V = H%. (7.111) 
Substituting (I7.110P into kinetic energy f)3.36p . we neglect the quadratic term 



Then 



j yoc 1*00 

Hi = - / VoVodu + / V 5i/jdu 



V (u,t). (7.112) 
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Here Vq(u, t) is a diven "slow" function of u and t. 
Now 



U = P-{^\ (7.H3) 
and equation (II. 17ft takes a closed form: 

lm{zz} = V . (7.114) 

In the absence of a potential energy Vo can be put constant. In this case 
equation (17. 114|) is exactly Laplace Growth Equation. In presence of potential 
energy V should be at least a function of time. In presence of gravity 



V = -g, 

V = —gt + constant. (7.115) 

This approximation to the Laplace Growth Equation derived-ed from hydrody- 
namics equation was obtained in Dyachenko &: Zakharov (1996)| , [Zakharov &: Dyachenko (1996)| . 



8 Cuts - possible type of solution 



Suppose both R and V have cuts along imaginary axis in the upper half-plane 

r b (t) r(s,t)ds rr/ , /•*>(*) v(s,t)ds 
R(w,t) = l- K ' ; , V(w,t)= v - ; (8.116) 

Ja{t) S + iw Ja(t) S + iw 



Here r(s,t) and v(s,t) are real functions of real argument s given in the 
interval [a(t),b(t)), and 

r(a, t) = r{b, t) = v(a, t) = v(b, t) = 0. (8.117) 

We show now that the anzats (18.1161) is the solution of the equations (13. 49 p . 
Let us calculate complex transport velocity U. 

U = P~ [RV + RV) . 
Anzats (I8.116P results in the following expression for U, 

L r_ j h v M dc f b f b r(s',t)v(s,t) + r(s,t)v(s',t) ds'ds 
s + iw J a J a s' + s s + iw 



20 



It is convenient to introduce the notation: 



G(s,t) = j'^^ds'. (8.119) 
Then U is very similar to V, namely 

U= f b < s ^ ds . (8.120) 

J a S + ilV 

Here 

u(s, t) = v(s, t) - F(s, t)v(s, t) - r(s, t)G(s, t) (8.120) 

It is obvious that u(s,t) is also equal to zero at the ends of the interval 
[a(t), b(t)}. Similar formula one get for B: 

B = P- \VV\ = t ^^-ds (8.120) 

I J J a S + IW V ' 

Here 

b(s,t) = v(s,t)G(s,t) (8.120) 

And again, b(s,t) is equal to zero at the ends of the interval [a(t), b(t)]. 

Now, let us calculate r.h.s of the equation for R, 

i (UH - U'R) = t f As,t)r(s\t)- r>(s,t)u(s>,t) ds , ds _ f'^M = 

J a J a [S + iW)(s' + IW) J a S + iw 

b r b u'(s, t)r(s', t) + u'(s', t)r(s, t) - r'(s, t)u(s ; , t) - r'(s', t)u(s, t) ds'ds 

J a S' — S S + iw 

If we introduce functions 

r ( S ') A, 



A(s,t) 


= P.V. [ 

J a 


C(s,t) 


= P.V. [ 

J a 


P(s,t) 


= -P.V. 


Q(s,t) 


= P.V. f 

J a 



v(s') 



r b r(s')G(s',t) + v(s')F(s',t) 

J a S 1 — S 

b v(s')G(s',t) 



s' — s 



ds (8.120) 



then equation for R reads: 



r(s, t) = -u'(s, t)r(s, t) - ii'(s, t)r(s, t) + r'(s, t)u(s, t) + r'(s, t)u(s, t) + u'(s, t) 

(8.120) 
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Doing the same thing for the equation for V one can obtain: 

i)(s, t) = u(s, t)v'(s, t) + u(s, t)v'(s, t) + b'(s, t)r(s, t) + b'(s, t)r(s, t) — b'(s, t) 

(8.120) 

or with new notation: 



r(s, t) + ((1 - A)G -C- P)r'(s, t) - (1 - A)(l - F)v'(s, t) = 

= ((A - 1)G' - A'G -C- P')r{s, t) + ((1 - F)A' - (1 - A)F')v{s,t), 

v{s, t) + ((1 - A)G -C- P)v\s, t) = 

= (Q' - GC')r(s,t) + ((1 - F)C - (1 - A)G')v(s, t) (8.121) 



9 Equations for spectral density on the cuts for implicit equations 



Here we will consider the hydrodynamics equations for y and ip given in the 
implicit form (11. 1 5 [) . Actually we will use instead of y and if) complex functions 
z and $, and the equations (11.151) can be rewritten for them: 

zz' - zz' = (9.121) 

$z' - &z + P~{^z' - &z) + &z' - $'! + P + {$z' - &z) = 0. (9.121) 

Let us make the following hypothesis. Suppose that both z and $ are analytic 
functions in the lower half-plane and both have the only singularity in the 
upper half-plane, namely a cut on the imaginary axis: 

X(t) < v < oo, 

X(t) is some unknown function of time. Let us denote the coordinate along 
the cut by s, so that 

w = u + iv = u + is. 

Now we introduce the spectral density on the cut for conformal mapping z 
and potential $: 



z 



1 f°° p(s') 



-ds\ 



2m Jx s' + iw 

$ = / yv ; ds. 9.122 

2ni Jx s' + iw 

Let us calculate z on the cut from right and from the left or, another words, 
for 

w = is ± e, e — > 0. 

Then 

z ± = —P.V. [°° p(g/) ds' T -pis). (9.122) 
2ni Jx s' — s ±ie 2 
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Then one can see that 

p(s) = z~ — z^ 
Let us introduce two more functions: 



and 



Then 



z[s = 



2ni 

Here is some useful notations: 

d d . d 

dw ids ds ' 

zt = iA t , z t z u ^A t A s . (9.123) 



Now let us consider spectral density for potential $: For values of $ on the 
right and left sides of the cut one can easily derive the following relations: 



$ - $ + = 2(/>(s), 

+ $+) = J-P.V. /°° (9.124) 
2 V ; 2tt A s'-s y j 

And similar to what we just did for z, let us use the notation: 

9 (s) = ±P.V. Holds', 
yv ; 2tt Jx s'-s ' 

B(s) = ^- raids'. (9.125) 
Then the following relations hold: 

i rW)j. 



Hs) = ^ Holds' = B(s), 
K J 2tt Jx s' + s y h 



—{$- _ $+) = rr 
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!($- + $+) =2(3). (9.126) 



Now we can easily rewrite the equation for i 



p t {l + A s ) - p s A t = -<j> s . (9.127) 

To write down the equation for $ (instead of fl9])) one should use the following 
useful formula: 

A-B- - A + B + = -[(A" - A + )(B~ + B + ) + (A~ + A + ){B- - B + )\. 
The the equation Qj takes the form: 



0t(l - fs) + ft<j>. + Pt9s ~ 9tPs + PtB s - Ps B t = 0, (9.128) 
or if we introduce new function 

C(*) = g{s) + B( S ) = -p.v. r ^K ds ' 

it is equal to: 

<M1 - fa) + Ms + PtC s - p s C t = 0. (9.129) 



These equations (I9.127P and (I9.129[) are the basic equations describing the 
evolution of the densities on the cut. 

Let us write them in the divergent form: 



j t [0(1 - /.) + P C S ] = §- s (-<f>ft + pC t ) (9.130) 

Assuming that spectral density of the potential and conformal mapping are 
equal to zero at the ends of the cut, we obtain the integral of motion (for mass 
or the fluid and its vertical momentum): 



/ p(l + A s )ds = M, 
/ [0(1 — f s ) + pC s ] ds = P. (9.131) 
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For the kinetic energy % 

i f°° - - 
U=-\ ($ + $)($' -$')dw, (9.131) 

8 J -oo 

one can easily obtain the following formula: 

1 r°° 

H = -j <j)A s ds. (9.132) 
10 Approximation for narrow cuts 



Suppose that cuts for R and V are far from the real axis, namely their width 
(b — a) is much less then the distance to the real axis 

{b — a) « a. 

Then one can approximate 

U = P-(VR + VR) 

as 

U ~ V C R + VR C - V c , (10.132) 
here V c is the value of V at the some point on the narrow cut of R, and R c is 
the value of R on the narrow cut of V. The last term in (jlOj) appears due to 
asymptotic of R at infinity. 

Here is the ground for this approximation. Both R and V have singularities in 
the lower half-plane, at the complex conjugate points with respect to R and 
V. If we consider narrow cut (at the same place for R and V), that means 
we assume V c± V c and R ~ R c being time dependent only. This assumption 
allows us to get the approximate expression for B also: 

B = P~(VV) ~ V C V. (10.132) 

It should be mentioned here that for the limiting case of infinitely narrow cur 
(it is nothing but pole) the approximation is exact. 

Substituting (jTOj) and (flO]) into the (13.491) we end up with the following equa- 
tions: 

R + iV c R' = iR c {VR' - V'R) (10.132) 
V + iV c V = iR c (VV). (10.132) 

In the moving framework 

t 



X = w — i / V c dt 
Jo 
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the equations (jTOj) and ( ITOj) read: 

R = iR c (VR' - V'R) (10.132) 

V = iR c (W), (10.132) 
where space derivative is now with respect to \- 

It is remarkable that we derive complex Hopf equation (fTOl . If we introduce 
here new time r(£), so that 

t(t) = R c (t) (10.132) 

and 

r V 

x = w -i / -^dr. (10.132) 

JO -R c 

Recall that R = than we finally get the following set of quadratic equations: 

z T = iVz (10.132) 

V T = iW. (10.132) 
These equations are Hamiltonian ones, with the Hamiltonian 

i r°° 

H = - V 2 z'd X (10.132) 
and canonical variables z and V, 

, = § , W = -f (10,32) 

Note that V and z' are analytic function in the lower half-plane, and TL is 
equal to zero. 

Equations ( ITU]) and (fill]) can be solved by the method of characteristics. 
Let us consider the following initial value problem for equations (110]) and (JTUJ: 



V(x,r)\ T=0 1 



^(x,r)| T=0 = X- (10.133) 

Here A and A are real positive constant. General solution of the Hopf's equa- 
tion (ITU]) is given by: 

F(V(x,t))=ix-tV( X ,t). 
Initial condition (110.1331) defines the function F(V))\ 

F(V) = £ - A, 
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and we end up with the quadratic equation for V(x, t): 

tV 2 - (A + %x)V + A = 0. 
Solution of this equation that satisfy the initial conditions is equal to: 



V( X , r) = Vl 2r (10.133) 

The branch of the square root in ffTOl) is chosen to provide zero asymptotic for 
V at infinity. 

General solution of (fTDj) with the velocity V(x, t) satisfying the Hopf 's equation 
is given by the formula 

z(x> r ) = G ( Z X - tV), 

with arbitrary function G. From the initial conditions f l 1 . 1 3 3 [) one can easily 
obtain that 

and for z(x, t) we get the expression: 



<X, r) = --{-A + i X + y/(X + ix) 2 ~ 4Ar}. (10.133) 
For R(x, t) one can get the following formula: 



2J(X + i X ) 2 -AAt 

R( X , r) = V , (10.133) 

A + %X + V( A + l X? ~ 4^4r 



Let us consider behavior of V(w,t) and R(w,t) just after the cut emerges, 
namely for 

A 2 

*«-. 

R(w,t) has two branch points 



iw\ = — A — / -f^-dr + 2v At, 
Jo R c 

i W2 = -\- / -JL(It-2VA^ (10.134) 
Jo R c 

Let us make an approximation to V c and R c . We will estimate value of V and 
R at the point which is complex conjugate to the branch point w\\ 
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1 A 

V c C.C. (V\ iw=—iw\) — 2"^' 
R c C.C. (R\iw=—iwi) — 1- 



Branch point hits the real axes when w\ becomes zero, at 



\2 

r = r ~ 0.34—. (10.133) 



In the vicinity of w\, R(w, r) behaves like this: 

2 



R(w, r) ~ -—rJi(w - Wl ) (10.133) 
[At)* 



(The latter must be deleted) 

Now we can estimate R c (t). According to (jTOl) and (flOl) it is equal to the value 
of the complex conjugate function R at the point w ~ wx- Or, it is equal to 
the value of R at the point w = —W\. Thus, 



2\/l- 



i? c (r) = -y= , (10.133) 

1 - ^ + ^1-2^ 



Recall that "true" time £ is related to r through the equation (ITUj) we can 
write down differential equation for time: 

dt=-—^. (10.133) 

After integration of (flOj) one can get 

5 1 

v2 A 1 / r. / A^\ 2 / O. / 2 



'-5^S^( 1 -^)-( 1 -^>- ™ 

For the small t and r relation ffTUj) gives 

8A 2 

Rescaled time r makes sense only for 

A 2 



t ~ r H tttT 2 . 



T - To 4A- 
At r = r "true" time tends to some finite value t ; 

_ 13 A 2 
^ to ~ 10 4A 



28 



and in the vicinity of To the following relation takes place: 

( To _ T )~iL(t -t) 2 . 

(The former must be deleted) 

Of course, the 'approximation of narrow cut' is valid only for the time 

T « T , 

but formally it can be extended to To, and one can obtain the limiting value 
of R(w, T ) from (TlOl) or (fT0|): 

R(w,t ) ~ JjViw. 



At the surface (w = u) 



X 

The slope of the surface is given by the ratio: 

dy Im{i?} 



u 

R(w,t ) ~ \ /V(l + *sign(u)). (10.132) 



dx Re{R} 



-sign^wj 



It means that the 90° angle appears at the surface. Again, this angle-type 
singularity takes place only for approximate equations (fit)]) and (fTOj) . 



11 Dynamics of zeroes and poles 



In this section we consider evolution of another initial condition for R, namely, 
instead of initially flat surface we will consider 



R(w, 0) = 1 + , (11.133) 

v ' b + iw y ' 

and the same (110.1331) for V. R(w, 0) has simple pole at w — ib and simple 
zero at w = i(a + b). Also it has right asymptotic behavior at infinity. Here 
we make important suggestion for the latter: pole of R is far from the cut 
of V It is obvious that this pole does not produce new singularity in V, and 
the approximations ffTO]) and ([TO]) hold ( at least for some time). The only 
difference is that for U additional term appears in ( JTUj) . namely: 

u^u+ const 
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12 Self-similar solutions 



Solution (fTUj) and ( TTUj) are self-similar for the variable: 
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A Poisson's Bracket 



Let us consider Poisson bracket between Hamiltonian "H and some functional 

a(z, z, $, $). 

r f { Set . 5a. 5a ■ 5a ^\ , . . „ s 

{ a,«} = /^-, + + + -4 jdu (A.0) 

Every part of the above integral we calculate separately. 

wH£}M'£}-*t£})* 
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Easy to see that {H,'H} = 0. 



B On the cubic equations 

C Hydrodynamics in the electric and magnetic field 

To be written. 

D Useful formulas 

y = Ex 
(x' + Ey'-)(x' - Hy'-)f — Jf 
E(A ■ B) = H(HA ■ EB) + AHB + BE A 

fj X ' M y ' 



P.V.r^ S, - a)ib - S,) ds' = n( a -±±-s) 

J a S'-S V 2 ' 

f b J(s'-a)(b-s>) a + b , 

/ — v— a ds = n—z — + s ~ y( b + s)(a + s)) 



A~B- - A + B + = -[(A- - A + )(B- + B + ) + (A~ + A + )(B- - B + )\. 
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